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affects  ecosystem  structure  as  drastically  as  the  extinction  of  a  bio¬ 
logical  species; (e.g. ,  Paine,  1966)1  As  demonstrated  by  public  concern 
about  projects  such  as  the  TVA  Tellico  Dam  Project  where  the  persistence 
of  the  snail  darter  Percina  tanas i  is  vascillating  (ftoTttenv  W7), 
extinction  of  a  biotic  component  of  a  food  chain  is  topical  as  well  as 
ecologically  fundamental. 

Viewing  the  importance  of  persistence-extinction  phenomena,  it  is 
somewhat  surprising  that  the  modelling  and  analysis  of  ecosystems  have 
not  concentrated  foremost  upon  species  survival.  The  Influence  of  stabil¬ 
ity  eventuation  in  the  engineering  and  mathematical  communities  was, 
undoubtedly.  Initially  too  Inviting  and  overwhelming. 

Ecosystem  stability  considerations  are  maturing  and  evolving 
with  practicality  and  realism  in  focus  (Holling,  1973;  May,  1974;  Maynard 
Smith,  1974,  Patten,  1974);  however,  a  general  applicable  theory  of 
stability  for  ecosystems  is,  at  best,  in  its  Infancy  stage. 

Goodman  (1975)  states  that  "minimally,  stability  means  persistence." 


In  this  article,  we  initiate  a  study  of  this  fundamental  aspect  of  eco¬ 
system  stability,  persistence-extinction,  with  emphasis  upon  development 
of  a  general  theory.  Persistence  attributes  of  simple  deterministic 
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food  chains  of  arbitrary  length  and  modelled  by  Lotka-Vol terra  dynamics 
are  determined.  Here  a  simple  food  chain  has  a  single  species  composing 
each  trophic  level  with  its  dynamics  governed  by  the  levels  immediately 
preceding  and  succeeding.  Persistence  is,  in  general,  a  global  property 
of  a  dynamical  system;  it  is  not  solution  structure  or  configuration 
oriented.  Equilibrium  analysis  used  In  conjunction  with  linearization 
techniques  has  been  a  principal  tool  used  for  studying  the  survival  problem 
(Freedman  and  Waltman,  1977;  May,  1974;  Resclgno  and  Jones,  1972).  Such 
approaches  have  not  been  particularly  fruitful  for  higher  dimensional 
food  chains  basically  due  to  Intrinsic  complex  analytical  and  topological 
problems;  for  example,  classification  of  recurrent  solutions  of  higher 
dimensional  dynamical  systems  is,  at  best,  difficult.  Other  techniques, 
for  example,  graph  theoretical  ones,  have  been  of  some  utility  In  higher 
dimensional  food  chains  (Yorke  and  Anderson,  1973). 

The  qualitative  approach  employed  here  introduces  the  concept  of 
a  persistence  or  extinction  function  which  Is  essentially  an  appropriate 
system  transformation  in  the  Liapunov  tradition.  This  technique  yields 
(global)  persistence  or  extinction  results  from  the  system  structure  with¬ 
out  a  priori  Information  of  the  asymptotic  character  of  the  model  solutions. 
The  conclusions  are  sharp  and  depend  only  upon  the  model  parameter  and 
food  chain  length.  , 

PRELIMINARY  TERMINOLOGY  AND  TECHNIQUES 

Let  R+  denote  the  nonnegative  real  numbers  and  r"  the  nonnegative 
cone  In  Rn;  that  is,  r"  «  ix  >  Rn  :  x  •  (x^,  ....  xn)  where  x^  f  R+, 

1  •  1.  2,  ....  n>.  The  positive  cone  In  Rn  is  R^  «  ix  *  R^  : 
x  •  (xr  ....  x„).  x,  '  OK 


# 
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The  models  considered  here  are  of  the  form 

Xj  ■  Xjf.j(Xj*  *  •  ••»  x^) i  1  ■  1 »  ». ,  ••••  n»  ■  d/dt;  (1) 

where  each  f^  Is  a  continuous  function  from  r"  to  R  and  Is  sufficiently 
smooth  to  guarantee  that  Initial  value  problems  for  equation  (1)  with 
the  Initial  position  x^(0)  '  0,  1  •  1,  2,  ....  n,  have  unique  solutions. 

The  primal  theme  develops  persistence  of  a  food  chain  modelled 
by  (1)  In  the  following  sense. 

Definition  1,  System  (1)  Is  persistent  If  each  solution  *  -  *(t) 

of  (1)  with  ^(0)  i  r"°  satisfies  11m  sup  *(t)  '  0  for  all 

t  *  t 

t  «  (0,T  1  where  [O.TJ  Is  the  maximal  Interval  of  existence  of  If 

(1)  Is  not  persistent,  then  there  Is  a  solution  with  >(0)  »  R™  where 

some  component,  say  satisfies  11m  *,(t)  •  0  for  some  i  In  (O.Tl. 

J  tM  J  * 

As  we  shall  be  concerned  with  persistence  of  the  complete  food  chain, 

any  solution  of  (1)  considered  subsequently  will  tacitly  have  Initial  value 

In  Rj°. 

If  It  Is  known  that  the  maximal  Interval  of  existence  of  all 
solutions  of  (1)  Is  [0,"'),  (that  Is,  T  •  >*)  then  It  follows  from  the 
assumed  uniqueness  of  solutions  of  Initial  value  problems  that  persistence 
of  (1)  Is  determined  by  using  only  v  •  »  In  Definition  1. 

In  the  sequel,  the  well  known  comparison  technique  Is  employed. 

It  Is  convenient  to  categorize  the  differential  equations  employed  as 
comparison  equations  In  the  following  terms.  Let  be  a  continuous 
function  from  R4  Into  R.  The  differential  equation 
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u'  ■  w(u)  (2) 

is  of  persistent  type  provided  any  solution  *  4>(t)  of  (2)  with 

4>(0)  >  0  satisfies  11m  sup  ^(t)  >  0  for  all  t  e  (0,®].  The  equation 

t  -*•  t 

(2)  is  of  extinction  type  provided  any  solution  v  *  v(t)  of  the  Initial 
value  problem  (2)  with  v(0)  >  0  satisfies  11m  v(t)  =  0  for  some  t  in 

t-»T 

(0,®].  For  example,  if  the  right  side  of  (2)  is  w(u)  =  au,  then  (2) 
is  of  persistent  type  whenever  a  ^  0  and  of  extinction  type  if  a  <  0. 

The  concept  of  a  persistence  (extinction)  function  is  now  intro¬ 
duced.  It  is  tacitly  required  that  the  functions  p  and  e  used  subsequently 
be  continuous  functions  from  r"  to  R+  which  are  continuously  differentiable 
on  r"°. 

Definition  2.  A  function  p  is  called  a  persistence  function  for 
system  (1)  if  the  following  are  satisfied: 

(i)  p(x-|,  x2,  ....  xn)-*0  If  xi  •*  0  for  some  1,  i  s  1,  2,  n; 
(1i)  p  satisfies  the  differential  inequality  p  >_w(p)  wherein 


’  xn^  5  ^  3x .  xlVxr  x2* 


(3) 


and  the  associated  comparison  differential  equation  u’  »  w(u)  is  of  per¬ 
sistent  type. 

Definition  3.  A  function  e  is  an  extinction  function  for  (1 )  if 
the  following  are  satisfied: 

(ill)  e(x<| ,  x2,  ...»  xn)  -*  0  only  if  some  x^  -*■  0,  1  ■  1,  2,  ....  n; 


i 
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(1 v)  e  satisfies  the  differential  Inequality  e  <_  u (c)  where  e 
Is  defined  as  in  (3)  and  the  associated  comparison  differential  equation 
u'  *  ui(u)  Is  of  extinction  type. 

The  existence  of  a  persistence  (extinction)  function  for  (1)  implies 
the  persistence  (nonpersistence;  l.e.,  extinction  of  a  component  species) 
of  system  (1). 

Theorem  1.  Let  p  be  a  persistence  function  for  system  (1).  Then, 

for  any  solution  4>  =  (^ ,  $2,  •••»  4>n)  of  (1)  with  maximal  Interval  of 

existence  [0,T  ),  lim  sup  <i> . (t)  >  0  for  each  t  e  (0,T  ]  and  each  i, 

9  t  ■+  T  9 

1  *  1,  2,  ....  n;  that  is,  system  (1)  is  persistent. 

Example  1.  The  function  p(x.|,x2)  s  x,x2/(l+x, )  is  a  persistence 
function  for  the  symbiotic  model 

X1  “  xl{!  "  IX1  +  2x2  +  xlx2  '  xl^’ 

*2  »  x2(l  +  x1  -  2x2); 

here,  p  satisfies  the  equation  p  *  p{|-  -  p)  which  is  of  persistent  type 
and  the  condition  p(x^,x2)  -*•  0  If  x^  -*■  0  follows  by  noting  that  x-j  ■+  0 
and  x2  -*>  «  are  incompatible. 

The  extinction  analogue  of  Theorem  1  is  our  next  result. 

Theorem  2.  Let  e  be  an  extinction  function  for  (1).  Then,  for 

any  solution  $  ■  U-j,  4>2,  ....  $n)  there  exists  an  1 ,  i  «  1 ,  2,  . . . ,  n 

and  a  t,  0  <  t  <  »,  such  that  11m  4>,(t)  *  0;  that  Is,  (1)  Is  not  persistent 

t+x  1 

and  extinction  of  some  component  species  results. 
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Proof  of  Theorem  1.  Suppose  that  system  (1)  Is  not  persistent;  then, 
there  exists  ate  (0,T  ]  and  a  j,  1  <  j  <  n,  such  that  11m  sup  *,(t)  <  0. 

*  ~  ~  t  -  T  0 

Condition  (1)  leads  to  the  conclusion  that  11m  p(*(t))  ■  0.  This 

t-*T 

conclusion  Is  not  consistent  with  the  assumption  that  the  comparison 
equation  u'  =  w(u)  has  no  solution  which  approaches  zero  since,  from 
elementary  differential  Inequalities,  it  follows  that  p($(t))  2iu(t), 
where  u(t)  Is  the  solution  of  the  initial  value  problem  u'  *  u>(u), 
u(0)  *  p(p(0)). 

The  proof  of  Theorem  2  is  similar  to  that  above  and  is  omitted. 
Example  2.  The  system 

x^  *  x^(l  +  ax-j)  a  •>  0 

x£  *  -x2(l  +  bx2)  b  >  0 

shows  that  considerations  of  the  maximal  interval  of  existence  in  the 
above  theorems  are  necessary.  The  function  p(x^,x2)  *  x^x2/(l  +  ax^  +  bx2) 
satisfies  p  *  0  and  Is  a  persistence  function.  Solutions  of 
x2  »  -x2  (1  ♦  bx2)  with  positive  initial  conditions  approach  zero  as 
t  approaches  Infinity;  however,  solutions  of  the  system  exist  only  on 
a  finite  Interval. 

As  will  be  demonstrated  below,  the  persistence  of  (1)  can  sometimes 
be  verified  even  If  the  persistence  function  does  not  satisfy  the  differen¬ 
tial  Inequality  throughout  the  entire  region  but  only  a  slab  of  the 
form  {(x^ ,  x2,  ....  xp)  e  r"°,  0  <■  Xj  ^  \jl  for  some  j,  j  ■  1 ,  2,  ....  n 


and  some  a.  >  0.  The  conclusion  in  these  circumstances  is  the  establish- 
ment  of  a  component  hierarchy  structured  by  the  dominance  of  the  species 
which  survive  the  struggle  for  existence. 

To  illustrate,  in  an  elementary  setting,  our  technique  for 
determining  species  survival,  we  consider  the  quadratic  model 

x-j  *  x.j(a  +  bx1  +  cx2) 

x£  =  x2(e  +  fx1  +  gx2) 

where  a,  b,  c,  e,  f,  and  g  are  constants  (not  necessarily  positive). 

Define  p(x, ,x2)  *  x-|X2;  then 

p  *  pC (a  +  c)  +  (b  +  f)x1  +  (c  +  g)x2]. 

Suppose  we  know,  a  priori,  that  solutions  are  bounded.  If  a  +  e  >  0, 
b  +  f  >  0  and  c  +  g  >  0  then  p  >_  (a  +  e)p  and  it  follows  that  the  system 
is  persistent.  Whenever  a  +  e  >  0  and  b  +  f  >  0  then  species  x2  survives 
since  p  >  0  on  the  slab  {(x^,x2)  :  x^  e  rJ®,  0  <  x2  <  (a  +  e)/(c  +  g)}. 
Similarly,  when  a+e>0andc+g>0  species  x^  survives.  These 
results  are,  admittedly,  not  sharp.  They  are  presented  only  to  demonstrate 
a  technique  utilized  subsequently  in  a  more  elaborate  setting. 

L0TKA-V01TERRA  FOOD  CHAINS 

By  employing  the  techniques  developed  in  the  previous  section, 
we  now  classify  persistence  in  a  simple  food  chain  modelled  by  Lotka-Vol terra 
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dynamics.  The  results,  which  depend  upon  the  interspecific  and  intra¬ 
specific  interaction  coefficients  and  the  length  of  the  predator-prey 
system,  determine  persistence-extinction  up  to  a  single  parameter 
(bifurcation) value.  The  presence  or  absence  of  a  carrying  capacity 
for  the  food  chain  resource  (lowest  trophic  level)  has  an  interesting 
effect  upon  the  character  of  the  persistence  criterion.  The  model 
considered  here  Is 


x'  “  x  (-a„n  +  a„  „  ,x„  ,), 
n  n'  no  n,n-l  n-1 


In  (4),  all  a^  are  positive  constants  with  the  exception  of  a^  which 
Is  nonnegative.  Before  proceeding  with  our  main  result  for  (4)  with  nonzero 
carrying  capacity  we  note  that  limitations  on  the  resource  level  are  pro¬ 
pagated  throughout  the  food  chain. 

Theorem  3.  All  solutions  of  (4)  with  positive  initial  conditions 
are  bounded  provided  a^  '0. 

Proof.  From  the  first  equation  in  (4),  It  follows  that 
xi  <  xl^alo  “  allxl^*  comPar*son  principle  establishes  the  bound 


x^(t)  <_  max  ix^O),  a1Q/a^}  for  t  c  R+. 
We  now  define  the  function  u  by 
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where 

n-1 

m  -  min  a  n,  b  »  max  1  x, (2a1n  -  anx,);  ]  j  a.  .  .  . 
lj^n  JU  1  10  n  '  k-1  K+u* 

Application  of  the  comparison  principle  leads  to 

u(t)  u(0)  exp  (-mt)  ♦  b/m; 

this  implies  x^,  i  *  2,  ....  n,  is  bounded  and  completes  the  proof  of 


the  theorem. 
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For  system  (4),  the  next  theorem  classifies  persistence  In  terms 
of  a  single  parameter  u.  Writing  the  food  chain  length  n  as  2m  +  1 
if  n  is  odd  or  2m  if  m  is  even,  we  define 


u  *  a 
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m 

y 

j=i 


*2 j  ,0 


a2i-2.2i-l 
isl  a2i ,21-1 


m 

V 

i=l 


a2i+l,0 


■pj  a2i-l  ,2i 
i  =  l  a2i+l  ,2i ' 


Theorem  4.  Suppose  the  resource  level  in  (4)  has  a  positive 
carrying  capacity  (a^  >  0).  Then,  the  food  chain  as  modelled  by  system 
(4)  is  persistent  provided  u  >  0;  it  is  not  persistent  if  u  *•  0. 

Proof.  First,  we  review  some  appropriate  properties  of  the 
dynamical  system  (4).  Bounding  hyperplanes  of  R™  of  the  form 
Hj  =  ((xp  X2,  x3,  ....  xn)  c  R^  :  x^.  =  0)  are  invariant  manifolds  for 
(4).  The  intersection  of  the  *>-limit  set  of  a  solution  (which 
is  bounded  by  Theorem  3)  with  such  a  hyperplane  is  therefore  a  compact 
invariant  set  which  contains  a  minimal  compact  invariant  set  (Nemytskii 
and  Stepanov,  1960).  A  classical  result  of  Birkhoff  implies  that  every 
trajectory  in  this  minimal  set  is  recurrent.  When  extinction  of  some 
constituent  of  the  food  chain  represented  by  (4)  occurs,  its  model  equiva¬ 
lent  is  a  trajectory  approaching  a  bounding  hyperplane.  This  extinction 
trajectory  converges  as  t  -*  »  (since,  as  previously  remarked,  finite 
extinction  time  is  not  possible)  to  a  set  containing  recurrent  solutions. 

Suppose  now  that  u  ^  0  and,  for  the  purpose  of  contradiction,  that 
system  (4)  has  a  solution  *  =  (^,  ^  ...,  ^n)  with  ^(0)  £  R^  which 


11 


satisfies  lim  $.(t)  ■  0  for  some  minimal  index  j.  This  assumption  can 
t—  J 

be  rephrased  as  H.  C\  si  f  i.  We  claim  that  any  recurrent  solution  in 
J  9 

this  intersection  lies  in  some  minimal  subsurface: 

Hk°  ■  <<V  x2 . *n>  :  N  '  0-  '  -  (  -  k  -  '•  Xj  *  0.  i  1  k>. 

kO 

That  is,  the  solution  is  in  H  but  not  in  any  subsurface  on  which  s  0 
for  some  i,  1  i  i  i  k  -  1.  The  trivial  solution  (0,  0,  ....  0)  cannot  be 
such  a  recurrent  solution  since  it  is  of  hyperbolic  type;  hence,  it  is 
not  in  »1  .  Now,  if  the  claim  is  false,  there  is  a  recurrent  solution 

't 

,(t)  *  ( v i ( t ) ,  ....  v  (t)  satisfying,  for  some  indices  c  and  m, 
v,(t)  =  0  and  v  (t)  /  0.  If  c  is  the  largest  index  with  (t)  :  0  the 

<  U\  v 

function  ,  satisfies  the  differential  equation 

VC+1  3  \c+l  ("av+l,0  VC+1  "  a;+l,v+2  vc+2)- 

Thus,  v  0  for  all  t,  a  fact  which  is  not  possible  for  a  recurrent 

solution.  The  minimality  of  is  ascertained. 

We  shall  now  show  that  }  cannot  have  its  assumed  asymptotic 
behavior.  This  is  accomplished  by  showing  that  for  some  r.  '  0, 
i  *  1,  2,  ....  n,  and  for  some  \.  N  0,  j  =  k,  ...,  n, 

J 

n  r. 

p  =  ]  [  x.  is  a  persistence  function  on  the  hyperrectangular  slab 
i*l  1 

s  *  <<V  x2 . *n)  >'  :  0  '  Xj  <  \Jt  j  1  k'. 

The  asymptotic  condition  required  for  o  to  be  a  persistence  function  follows 
readily  from  Theorem  3  and  the  definition  of  .  It  remains  to  show  that 
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Is  a  solution  of  a  differential  Inequality  whose  comparison  equation  is  of 
persistent  type. 

From  (4),  we  obtain 


p  *  Crl(a10  "  allxl  "  a12x2)  +  X  rj('ajO  +  aj,j-lxj-l  '  aj »j+lXj+l 5 


+  r  (-a  «  +  a  ,  x  ,11 
nv  nO  n,n-l  n-l;j 


*  ptrla0  -  L  rjajO  +  (r2a21  -  rlall)xl 

+  X(rjaj,j-1  "  rj-2aj-2,j-l)xj-l  '  rn-lan-l  ,nxn^' 


Now,  choose  the  r . 's,  1  £  j  <_  k,  so  that 

J 


r2/ri  =  an/a2i  and 


rj/rj-2  =  aj-2,j-l/aj,j-r  j  =  3’  4 . k- 


Using  these  relationships  to  determine  the  quotient  r./r,  as  a  function 

J  * 

of  the  coefficients,  we  find  if  j  is  even  and  4  £  j  <_  k. 


*  !ll  ¥$  a2i-2,2i-l . 
rl  a21  1=2  a2i ,2i-l 


while  if  j  is  odd  and  3  <_  j  k, 


Substitution  of  (6)  Into  (5)  leads  to 


p  «  •J[r1wk  -  ^  rj*jO  +  (rjaj,j-1  ’  rj-2aj-2,j-l)xj-l 


-  rn  ,a  .  x  ] 
n-  I  n- 1 ,n  n 


(7) 


where  If  k  Is  odd  (k  *  2p  +  1)  or  even  (k  *  2p),  Is  given  by 

?  ,  A  a21-2.21-l 
wk  *10  “  ^  *21  0  I  I  * 

k  iu  jb1  41, u  ia}  *21,21-1 


l 


*»1  2?+1’°  1-1  *21+1 ,21 

Since,  by  hypothesis,  u  Is  positive,  uk  Is  positive.  To  establish 
that  p  >  0  In  S,  we  select  the  remaining  r^.l-k  +  l.k  +  2,  ...,n, 
and  the  Aj,  1  «  k,  k  +  1 ,  . . . ,  n  sufficiently  small  so  that  the  magnitude 
of  the  remaining  terms  Is  less  than  r^.  The  expression  in  the  square 
brackets  In  (7)  Is  positive  on  S  so  p  >  0  and  as  previously  remarked, 
extinction  cannot  occur.  This  contradiction  establishes  that  u  >  0  Is 
sufficient  for  the  persistence  of  (4). 

To  demonstrate  that  u  «•  0  Implies  extinction,  the  extinction  function 
n  r, 

c  ■  e(x,  ,  x0,  ....  x  J  ■  TT  x4  .where  all  of  the  r.'s  are  chosen  according 
1  ‘  n  1  1 

to  the  scheme  (6)  used  above, Is  employed.  Then,  from  (4),  It  follows  that 


r 


t  satisfies  the  inequality 


G 


£tri"  ■  r.-iVi,.,»!i?c' 


with  a  comparison  equation  of  extinction  type.  Theorem  2  implies  that 
the  food  chain  is  not  persistent. 

We  now  consider  a  food  chain  modelled  by  (4)  in  which  the  resource 

level  has  no  carrying  capacity,  that  is,  a^  *  0.  In  this  setting,  the 

persistence  parameter  Ug  depends  more  intricately  on  the  dimension  of 

the  system  than  the  previous  case  with  nonzero  carrying  capacity.  If 

the  food  chain  length  n  is  either  2m  +  1  or  2m  +  2  for  some  positive 

m  j  a^j  i 

integer  m  then  uQ  -  a]Q  -  [  a2i+1  Q  IT 

u  iu  j=1  a2i+lj2i 


Theorem  5.  Let  a^  a  0.  The  food  chain  modelled  by  system  (4) 
is  persistent  provided  uQ  >  0;  it  is  not  persistent  if  Ug  <  0. 

Proof.  Theorem  3,  which  established  the  boundedness  of  solutions 
of  (4)  assuming  a^  '0,  was  required  in  the  proof  of  Theorem  4.  The 
analogous  lemma  employed  here  is  that  any  solution  *  of  (4)  with 
*(0)  e  such  that  lim^  ^ ^ ( t )  ■  0  for  some  j,  is  bounded.  To  establish 
this,  we  first  note  that  1  (t)  *  0  implies  1  ^ *  0  for 

1  i  j.  From  (4)  it  follows  that  for  t  sufficiently  large 


(t)  - 


Thus,  11m„  * 

r-wn 


4  +  1 


*j+l(t)t’\j+l,0  +  Vl^j^  "  aj+l,j+2'l'j->-2(t)] 

-  ♦j+i(t,c‘aJ+i.o/2]- 

(t)  ■  0;  an  induction  completes  the  argument. 
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It,  then,  suffices  to  prove  the  assertion  in  the  case  that  4> 

is  such  that  1 imt^,^n( t )  ■  0.  To  this  end,  we  use  the  classical  Lotka- 

n 

Volterra  auxiliary  function  v{xj,  x2,  ....  xn)  ■  l  (x^-e^-i^ in  x^/6^), 
where  o^,  ^  are  positive  parameters  to  be  chosen  later.  Boundedness  of 
$  will  follow  if  it  is  shown  that  v(t)  ■  v($(t))  is  eventually  non increasing. 
Employing  (4)  with  some  rearrangement  leads  to 


n-1 


v(x.j,  x2,  ....  xn)  •  ^  C'aiai ,i+l  +  a1+lai+l,1^xr!jl^xi+r‘5i+l^ 


+  VvM^lO  -  ai2^2^ 

n-1 

+  a^Xi  ’  *^'*10  +  *i  ,i-lBi-l  “  ai,i+l’l5i+l) 


+  "  O (~®„n  +  ®n  n  1  1^’ 

n  n  n  nu  n,n-i  n-1 


The  parameters  i  »  1,  ....  n  are  selected  so  that  the  first 
sum  of  n  -  1  terms  vanishes.  If  n  is  even,  it  is  possible  to  choose  the 
flj  so  that  each  of  the  remaining  terms  is  zero;  thus.  In  this  case, 
aj.iJj  exist  so  that  v  h  0  from  which  the  boundedness  of  $  follows. 
Whenever  n  Is  odd,  we  can  choose  the  6^  parameters  in  such  a  way  that  all 
terms  vanish  with  the  exception  of  the  last  one;  hence, 


■A  • 


6n)('*nO 


*  Vn-1  Vl 1 


aA  -  8nll,0 


A  *21*1.21. 

1-1  *21-1,21 
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Thfs  leads  to 


»<t)  •  Vo  ^  IT 1±Lli  <♦„<*>  -  6n> 
n0  1-1  a21 -1 ,21  n  n 


which,  since  Ug  >  0,  satisfies  v(t)  <  0  for  all  t  sufficiently  large. 
Therefore,  any  solution  which  does  not  persist  is  bounded. 

The  remainder  of  the  proof  is  the  same  as  that  of  Theorem  4  up 
through  the  calculation  of  p;  from  (5)  with  a^  *  0  we  obtain 


5  >  ptr,(«)0  - 

*  jj  ‘Vi.l-l  '  rl-2al-2,l-l)xl-l  ‘ 


(8) 


I 

Now  let  k  be  written  ask*fc  =  2p  +  lork*s.  +  l=2p  +  2 
where  i  is  an  odd  integer.  We  choose  J  *  1*  2»  . ••»  P  as  in  (6b) 


rA 


a2i-l,2i 

a21+l ,21 ' 


Since  >  0,  it  follows  that 


-  ?  a  ^21±1 

10  ^  a2j+l  ,0  r, 


0. 


Select  the  remaining  rgj,  j  *  1,  ...»  p  (or  p  +  1  if  k  ■  2p  +  2)  so  that 

k 

*10  ■  j£2  Vi  aj0  >  (9) 
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If  q  Is  the  largest  even  Index,  q  ^  k,  we  now  adjust  r2,  r^,  ....  rq 
sequentially  (to  bound  successive  terms  on  the  right  side  of  (8)  while 
preserving  Inequality  (9))  according  to  the  following  scheme.  Choose 
rq-2  P°ssibly  smaller  than  the  original  choice  so  that 


VV2  -  aq-2,q-l/aq,q-T 

In  a  similar  manner,  choose  rq_^  smaller  if  necessary  so  that 

rq-2/rq-4  -  aq-4,q-3/aq-2,q-3- 

Continue  adjusting  r2i  until  r2  Is  revised. 

We  now  have 

k  n 

5  i»Crl(*lo  '  ,l2  <ri/ri)aio>  '  riato 

*  1«k+l  ’  r1-2a<-Z.1-''Xi-l  "  VlVl.A1, 

Proceeding  as  In  Theorem  4,  we  can  select  the  remaining  r^'s  and  the  a^'s 
so  that  the  expression  In  the  brackets  remains  nonnegative  on  S.  This 
leads  to  the  desired  conclusion  that  uQ  >  0  Implies  system  persistence. 

The  proof  that  Pg  <  0  Implies  nonpersistence  Is  analogous  to  the 
similar  portion  of  the  proof  of  Theorem  4.  As  an  extinction  function  use 
c(x1$  x2,  ....  xn)  ■  x^^1  where  r^  are  determined  by  (6b).  Then, 

when  n  Is  odd,  the  differential  equation 


J 
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£  *  rluOc 

results;  If  n  Is  even,  we  obtain 

1  ■  £trl“0  -  rn-lan-l ,nxn^* 

from  which  nonpersistence  follows. 

COMMENT  AND  COMPARISON 

While  the  above  results  are  Independent  of  system  equilibrium 
considerations,  the  conclusions  have  a  positive  correlation  with  the 
location  of  the  equilibrium  points  when  they  exist.  For  convenience  in 
demonstrating  these  and  other  relationships,  we  restrict  consideration  to 
food  chains  of  length  four  with  nonzero  resource  carrying  capacity: 

xi  3  xi(aio  -  anxi  ■  ai2x2) 

x2  “  x2^"a2Q  +  a21xl  "  a23x3  ^ 

(10) 

x3  "  X3('a30  +  a32x2  '  a34x4} 
x4  *  x4^’a40  +  a43x3  ^  * 


The  equilibrium  points  of  (10)are  Indicated  In  Figure  1. 


Figure  1 


If  we  designate  the  1—  coordinate  of  equilibrium  point  I  by  xj,  II  by 
xj*,  etc.,  the  persistence  Theorem  4  may  be  phrased  as  (10)  is  persistent 
if  the  coordinate  equilibriums  satisfy  the  inequalities 


x*  <  xjV  (with  the  sequential  ordering  x][  <  xjV  >  xj11  <  xj*  >  0); 


or  equivalently. 


V  iv  V  IV  III 

x2  >  x2  (with  the  sequential  ordering  x2  >  x£  <  x2  >  0) 


or  equivalently. 


4  <  x'\  <  x'V  >  0); 


or  equivalently. 


V 

x4  >  0. 


The  last  inequality  is  just  a  restatement  of  the  persistence  inequality  y  >  0. 
Analogously,  if  the  coordinate  ordering  is  reversed,  nonpersistence  results. 

This  equilibrium  coordinate  ordering  indicates  that  Introduction 
of  a  new  top  predator  into  an  ecosystem  would  result  in  an  order  interchange 
of  equilibrium  populations  on  an  alternating  scheme.  For  example,  in 
an  aquatic  system  with  trophic  levels  composed  of  phytoplankton,  heri- 
vores,  and  primary  carnivores,  the  model  conclusion  predicts  that 
the  addition  of  a  secondary  carnivore  reduces  the  equilibrium  levels 
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of  the  phytoplankton  and  primary  carnivores  while  increasing  the  herbivore 
population. 

The  component  ordering  also  indicates  that  the  persistence-extinction 
mechanism  is  mathematically  one  of  bifurcation  associated  with  the  parameter 
u;  essentially  the  equilibrium  V  for  (10)  "splits"  off  from  IV  as  the 
parameters  transact  from  the  extinction  (u  <  0)  to  persistence  range 
(u  >  0). 

In  the  case  of  a  persistent  system,  the  resource  transfer  is  modu¬ 
lated  by  trophic  level  Influence  factors  which  determine  to  some  extent  the 
magnitude  of  persistence.  These  factors  are  indigenous  to  the  specific 
trophic  level  and  relate  persistence  in  the  food  chain  hierarchy.  If  the 
model  is  interpreted  as  nutrient  flow,  energy  flow,  etc.,  componentwise 
throughout  the  system,  trophic  level  1  contribution  factor,  a^/a^,  is 
the  persistence  multiplicative  effect  of  level  1  on  trophic  level  2. 

Trophic  level  2  has  a  multiplicative  contribution  factor  of 


■  ■  on  trophic  level  3  while  level  has  multiplicative  contribution 
alla23 


factor 


alla23a32 

a12a21a34 


for  level  4. 


The  antipodal  character  of  certain  of  the 


trophic  level  influence  factors  might  not  have  been  anticipated. 

In  general,  the  bifurcation  value  is  persistence  indeterminant. 

In  the  noncarrying  capacity  model,  Freedman  and  Waltman  (1977)  have  shown, 
for  the  three  dimensional  system,  u  ■  0  implies  persistence.  On  the  other 
hand,  for  the  four  dimensional  food  chain  (10)  with  a^  *  0  there  is  a  con¬ 
tinuum  of  critical  points  in  the  positive  x-^Xj  space  when  uQ  «  0.  A 
subset  of  these  satisfies  x3  <  a^/a^.  Linearization  about  any  of  these 
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latter  equilibrium  points  indicates  that  each  has  an  asymptotically  stable 
manifold  which  Intersects  the  positive  cone  in  (the  eigenvalue  in  the 
direction  is  negative).  Thus,  there  are  solutions  with  positive  initial 
conditions  which  tend  toward  these  equilibriums  as  t  approaches  infinity 
and  u  ■  0  corresponds  to  nonpersistence. 

Goh  (1977)  has  recently  shown,  under  certain  conditions,  that  when¬ 
ever  systems  of  form  (1)  have  an  equilibrium  in  the  positive  cone  in  Rn, 
it  is  globally  asymptotically  stable.  This  is,  of  course,  a  persistence 
result  when  it  it  applicable.  To  contrast  our  conclusions  with  Goh's, 
we  have  found  that  the  existence  of  an  equilibrium  in  the  positive  cone 
is  equivalent  to  persistence  modulo  the  bifurcation  value  u  *  0. 

It  Is  worth  noting  that  for  odd  dimensional  models  (4)  without  car¬ 
rying  capacity,  persistence  (uQ  >  0)  results  even  though  there  is  no 
positive  equilibrium. 

SUMMARY 

Persistence-extinction  in  simple  food  chains  modelled  by  Lotka- 
Vol terra  dynamics  is  governed  by  a  single  parameter  which  depends  upon 
the  interspecific  interaction  coefficients,  the  Intraspecific  interaction 
coefficients,  and  the  length  of  the  food  chain.  In  persistent  systems 
with  nonzero  carrying  capacity,  two  new  features  predominate.  Trophic 
level  Influence  factors  relate  persistence  on  different  trophic  levels 
and  determine,  In  conjunction  with  the  persistence  parameter,  the 
magnitude  of  persistence.  Equilibrium  component  ordering  which  results 
in  persistent  systems  mandates  once  again  that  systems  need  to  be  studied 


I 


i 


on  the  complete  ecosystem  level;  static  field  measurements  reflect  species 
location  in  the  food  chain,  the  total  length  of  the  food  chain  and  assume 
characteristics  according  these  factors. 
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LEGEND  FOR  FIGURE  1 

Equilibrium  points  of  (10)  are  listed  hierarchically  In  terms  of 
coordinates.  " ~  ^ "  represents  implications  for  the  appropriate 
coordinate  orderings.  Trophic  level  influence  factors  for  each  level 
are  listed  (see  text  for  explanation). 


